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Abstract 



The most important part of the order a] corrections to the charm compo- 
nent of the charged-current structure functions ^2(2;, Q'^) and ^3(0;, Q"^) have 
r^' been calculated. This calculation is based on the asymptotic form of the 

^ . heavy-quark coefficient functions corresponding to the higher order correc- 

tions to the W-boson-gluon fusion process. These coefficient functions which 
'k>( \ are in principle only valid for Q^ ^ m? can be also used to estimate the 

^ \ order a1 contributions at lower Q^ values provided x < 0.1. It turns out that 

the above corrections are appreciable in the large Q^-region and they explain 
the discrepancy found for the structure functions between the fixed-flavour 
scheme (FFS) and the variable-flavour-number scheme (VFNS). These cor- 
rections also hamper the extraction of the strange-quark density from the 
data obtained for the charged-current and the electromagnetic-current pro- 
cesses. 
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1 Introduction 

The study of charm production in deep-inelastic lepton-hadron scatterings 
provides us with important information about heavy-quark production mech- 
anisms. The latter enables us to extract parton densities from the deep- 
inelastic data in certain kinematical regions which are difficult to explore in 
other hard processes. The best example is the gluon density G(a;, /x^) where 
/i denotes the factorization scale. It can be measured at small x in deep- 
inelastic electroproduction e^ -|- A^ ^ e^ +'X' via the photon-gluon fusion 
process 

j* + g^c + c (1.1) 

where 'X' denotes any inclusive final state. The above process dominates 
electroproduction of the charm quark if one assumes that the probability to 
find a charm quark inside the proton is zero so that the fiavour-excitation 
process 

7* + c^c (1.2) 



does not occur. The reactions (1^) and ( |1.2|) are sometimes also referred 
to as extrinsic- and intrinsic-charm production. In the subsequent part of 
this paper we will assume that there is no charm quark in the initial state. 
Another parton density which can be measured via charm production is the 
strange-quark density s(x,/i^). The latter shows up in the charged-current 
process z/(z/) + p — > l~{l~^) +'X' (/ = e,/x) which in lowest order is given by 
the fiavour-excitation mechanism 

W^ + s^c, W' + s^c. (1.3) 

In next-to-leading order (NLO) we have the following parton subprocesses. 
The first one is represented by the gluon bremsstrahlung process 

W+ + s ^c + g, W~ + s -^c + g (1.4) 



which includes the one- loop corrections to the Born reaction (|1.3| ). Further- 
more we have the W-boson-gluon fusion mechanism given by 

W^ + g ^ c + s, W~ +g ^c + s. (1.5) 



Although for x > 0.1 charm production is determined by flavour excita- 
tion it turns out that the W-boson-gluon fusion process dominates the order 
Us corrections for x < 0.1 and constitues a large background for reaction 
( |1.3| ). This means that in the region x < 0.1 an accurate determination 
of s(a;,/i^) requires a good knowledge of G(x, /i^). Charm production in 
charged- current interactions was studied in the experiments carried out by 
the CCFR-coUaboration in 0], 0]. Besides the strange-quark density one 
has also determined one of the elements of the CKM (Cabibbo-Kobayashi- 
Moskawa) matrix, called \Vcd\i which is approximately equal to | sin^d where 
dc denotes the Cabbibo angle. However there exists an alternative way to 
determine the strange-quark density s(x,yU^). Here the latter follows from 
the relation 



A(x,/i2) = ^F2^^(x,/i2) -3Fr(x,/x2). (1.6) 

o 

Here F2^ and F2 are the structure functions corresponding to charged- 
current and electromagnetic-current induced processes respectively, where N 
is represented by an isoscalar nuclear target. Neglecting nuclear corrections 
one obtains 



Ff (x,/i2) = ^(^F^V,/^') +^2"(^,/i')j (1-7) 

with / = e^,/i^, z/g, v^ and p and n stand for the proton and the neutron 
respectively. Further F!^^ denotes the full structure function where besides 
charm-quark production also the other light-parton (u,d,s,g) subprocesses 
are involved. If the structure functions appearing on the right-hand side of 



L6D are computed on the Born level we obtain in a three-flavour scheme 



A(x, /i^) = a; s(x, /i^) . (1.8) 

In the case of a four-flavour scheme where now also the charm quark is 
represented by a parton density Eq. (|1.6| ) becomes 



A(x,/i^) = X s{x,jj,'^) — X c{x, fx^) . (1.9) 



The three- and four-flavour scheme are in the hterature very often referred to 
as fixed-flavour scheme (FFS) and variable-flavour-number scheme (VFNS) 
respectively. The structure functions F2^ and F2 have been measured by 
the CCFR collaboration ||^, [Q and the NMC-coUaboration ^. If one takes 
the data from 0, then the result obtained for s(x, /i^) is quite close to 
the one measured in charm production ^, p (see also the discussion in 0). 
However there exists a discrepancy between the data of @] and |Q which is 
not clarified yet. If we adopt the three- flavour scheme for the subsequent 
part of our paper the identity in ( [L.8|) only holds if the charm-quark mass is 
zero and the QCD corrections are neglected. However as we will see later on 
these corrections will considerably alter Eq. ( |1.8|) in particular in the small 
x-region. This is mainly due to the W-boson-gluon fusion process in ( |1.5| ) 
and its higher order QCD corrections. The corrections due to the mass of the 
charm quark become smaller when the virtuality of the W-boson denoted by 
Q^ will become very large. The order ag corrections given by the reactions 
( |1.4|) , ( |l75|) were calculated in (see also 0, f^). The dependence on the 
mass of the charm quark was studied in [Q and a comparison between the 
FFS [0 and VFNS [|ll| was made in [l^, |jl3|. In the latter references the 



dependence of the QCD corrected structure functions on the factorization 
scale was studied too. One of the main conclusions of these investigations is 
that the QCD corrections for x > 0.1 are small whereas for x < 0.1 the bulk 
of the corrections is constituted by the W-boson-gluon fusion process ( |1.5|) . 
These corrections become even larger when Q^ increases. In particular at 
HERA 0, [|15l, H where 200 < Q^ < 10^ (GeV/c)2, the size of the con- 



tribution of the latter process is responsible for the discrepancy between the 
FFS and VFNS in describing the charm component of the structure functions 
denoted by Fk^c {k = 2, 3) as discussed in [0. Therefore we will concentrate 
ourselves in this paper on the calculation of the order a^ contributions to 
Fkc which are due to the order a^ corrections to the W-boson-gluon fusion 
process in ( |1.5| ).This also includes the contributions from other heavy-quark 
subprocesses initiated by the light quarks u,d and s. Unfortunately for some 
of them we only have the asymptotic expressions of the heavy-quark coeffi- 
cient functions which are valid for Q^ ^ m'^ where m stands for the heavy- 
quark mass. This implies that strictly speaking they are only applicable to 
the HERA experiments where the characteristic values for Q^ are very large. 
However as is shown for the electromagnetic-current process [jl^-pOl , where 



the exact order a^ corrections to the photon-gluon fusion process are known 



np, the asymptotic heavy-quark coefficient functions can be used for much 
lower Q^ values measured at fixed-target experiments. If Q^ is an order of 
magnitude larger than vn? one gets a reasonable agreement within 10% be- 
tween the predictions obtained from the exact and asymptotic heavy-quark 
coefficient functions (see |ll8|-||20[|). 

Our paper will be organized as follows. In Section |^ we present the 
formulae for the charged-current deep-inelastic structure functions expressed 
in convolutions of the parton densities and the coefficient functions. Further 
we construct the asymptotic form of the order a^ contributions to the heavy- 
quark coefficient functions. In Section ^ we discuss our results and show the 
effect of the latter on the charm component of the charged-current structure 
functions. We also show how the extraction of the strange density from 
relation ( |1.6| ) is influenced by the above corrections. The derivation of the 
asymptotic heavy-quark coefficient functions will be presented in Appendix 
A. 



2 0(al) contributions to heavy-flavour coeffi- 
cient functions in charged- current interac- 
tions 

In this section we will present the asymptotic forms of the heavy-quark coef- 
ficient functions up tp order a^ which contribute to the structure functions 
F2{x,Q'^) and F3{x,Q'^) measured in charged-current deep-inelastic lepton- 
hadron scattering. Although these coefficient functions are only valid at 
Q^ ^ m'^, characteristic of the values of Q^ observed at HERA, they can 
be also used to make a reasonable estimate of the charm component of the 
structure functions measured in fixed-target experiments where Q^ is much 
smaller. Here and in the next section we will compute all structure func- 
tions in the so called fixed-flavour scheme (FFS). In the case of the charm 
production this means that the constituents of the hadron are only given by 
the gluon and the three light flavours u,d,s which are described by parton 
densities. The charm quark only appears in the flnal state of the parton 
subprocesses and its contribution to the structure functions is described by 
the heavy-quark coefficient functions. As is shown in ||2^ the FFS is suit- 
able to describe the charm component to the structure functions at small 
Q^-values. However at large Q^ where the large logarithms of the type 
ln*((5^/m^) ln-'(/i^/m^) dominate the heavy-quark coefficient functions this 
description is not adequate anymore. In this regime it is much more useful 
to adopt the variable-flavour-number scheme (VFNS) where now also the 
charm-quark contribution is described by a parton density. This parton den- 
sity represents the resummation of the large logarithms above in all orders 
of perturbation theory. We will comment on the VFNS approach in the case 
of charm production at HERA collider in the next section. In the FFS (here 
three-flavour scheme) the light-parton contribution to the structure functions 
measured in the process ui + p -^ I +'X' (/ = e, /x) is given by 

Fr^(x,Q2)^Fr(x,g2) = 



Jx Z 



^4prV?? 4 H-G f,.^ cii4 



■^ 2 1,7, „„„2//i N j/ -^ ,.2 



+ Ufc,«M( -,/i 1 +&A:,dCOS (6'c)c?l -,yU' 




(2.1) 



a2{x) = 2x, asix) = 2 (2.2) 

and for i = u,d, s; i = u,d,s 

b2,i = hi = 1, hn = -h'i = 1 • (2-3) 

Here G{z, fi^) denotes the gluon density and the singlet combination of parton 
densities in the three-flavour scheme is given by 

S3(z,yU^) = u{z, f/) + u{z, fi'^) + d{z, ^'^) + d{z, fi'^) 

+s{z,fx'^) + s{z,f/). (2.4) 

Like the parton densities the coefficient functions Ck^i {k = 2,3;i = q, g) can 
be distinguished into a singlet and non-singlet part indicated by the super- 
scripts S and NS respectively. In ( p.l|) we have split the singlet coefficient 
function Cf in the following way 

^k,q — ^k,q + ^k,q ■ (2-5) 

The purely-singlet part of the coefficient function indicated by PS originates 
from quark subprocesses where the projection on the non-singlet channel 
yields zero, so that only singlet contributions remain. They are character- 
ized by those Feynman graphs in which only gluons are exchanged in the 
t-channel. Such graphs show up for the first time in order a^. Furthermore 
because of charge conjugation invariance of the strong interactions one has 
the property 

Cl^ = 0, Cl = , (2.6) 

so that -F3 is determined by the non-singlet combination of parton densities 
and Cfg only. In the case of F2 we have extracted the overall dependence 
of the coefficient functions Cf"^ and C2g on the number of light flavours nf 
so that we will denote them by C|[g and Cf^^ respectively. Finally, in the 
three-flavour scheme, the structure functions depend on the Cabbibo angle 
9c only. 

7 



The charm contribution to the structure functions F^^ {k = 2, 3) can be 
written as 



' X 







+bk,sCos\9c)si-,f/] \H, 



(2.7) 



with 



Zth 



Q' 



Q2 j^ 4^2 • 



^th 



Q' 



Q2 _j_ ^ 



(2.8) 



The above upper boundaries of the integrals in ( p.7|) follow from the threshold 
conditions. In the first part of Eq. ( |2.7| ) the heavy-quark coefficient functions 



Lj:^ originate from parton subprocesses with a heavy quark and a heavy anti- 
quark in the final state. Here the W-boson is attached to the light quarks 
u,d,s only. Hence the partonic centre of mass energy squared s has to satisfy 
the condition s > Am? with s = Q^(l — z)/z. In the second part of (|2.7| ) the 



heavy-quark coefficient functions if^ originate from the parton subprocesses 
where either the heavy quark or the heavy anti-quark appears in the final 
state which is accompanied by a light (anti-) quark. Furthermore the W- 
boson also couples to the heavy (anti-) quark. Therefore s > iv?. Notice 
that in ( |2.71 ) we have limited ourselves to those parton subprocesses which 
have no more than two heavy (anti-) quarks in the final state. The notation 
in ( p.7|) is the same as given for the light-parton contributions to the structure 
functions listed below ( |2.1| ). Because of charge conjugation invariance of the 



strong interactions the L^^ satisfy similar relations as quoted in (p.6|) i.e. 



^rr=o, L^/=o. (2.9) 

The structure functions F^^^ and F^^ which show up in the anti-neutrino 
reaction ui + p —>■ I +'X' can be derived from Eqs. ( |2.1| ) and ( |2.7| ) by the 
replacements u ^->- u, d ^->- d, s ■^^ s. Moreover in the case of F^^ the heavy- 
quark coefficient functions H^^ and H^g get a relative minus sign with 



respect to the ones showing up in F^^ (|2.7|). 



In the case of an isoscalar target one also has to compute the structure 
functions for the neutron given by F^"", F^", F^", F^". The latter can be 
derived from the proton structure functions by the replacements u ^^ d and 
u ^^ d. 

In the next section we want to compute the structure functions in (|2.1| ) 
and ( p. 71) up to next-to- leading order (NLO). Moreover we want to include the 
order a^ contributions to the heavy-quark coefficient functions Lf^'^ , Hj^j^ 
and Hj^g appearing in F^^^ (^.T]). Notice that up to order a^ the coefficient 
functions L2^q and L2^g do not contribute because the corresponding parton 
subprocesses are of the order a^ and a^ respectively. For L^^ and L^g see 
( |2.9| ). Furthermore we do not include the order af contributions to if^^ . 
The latter coefficient function is determined by the higher order QCD cor- 
rections to the flavour-excitation process ( |1.3| ). In order as it turns out that 
these corrections are much smaller than those originating from the W-boson- 
gluon fusion process ( P!75| ) leading to the coefficient function H/^g . Hence 
we expect that beyond order as the radiative corrections are dominated by 
the W-boson-gluon fusion mechanism which is also indirectly present in the 
parton subprocess contributing to Hf^^ . 

For the computation of F^^, F^^ ( |2.1| ) one needs the next-to-leading-log 
(NLL) parton densities and the light-parton coefficient functions Ck^i {k = 
2,3;i = q,g) corrected up to order as- Notice that C|^^ vanishes up to this 
order. The representation in the MS-scheme can be found in Appendix I of 



2^. For the calculation of F^^, F^^ the heavy-quark coefficient functions 



Hj^g and Hj^ g have been computed in up to order a^ (see also p, |§|). 
The latter are represented in the MS-scheme in Appendix A of [^. In this 
paper we want to compute the order a^ contributions to Lj^'^ , if^ ^ and 
H}^ g . The coefficient function Lj, ^ is given by the Compton reaction 



W + q{q)^q'{q') + c + c (2.10) 

where the heavy-quark pair originates from gluon sphtting and the gluon is 
radiated off the hght quarks q and q' (g, q' = u, d, s). The reaction above also 
includes the two-loop vertex correction to the subprocess W + q{q) — » q'{q') 
containing the charm-loop contribution to the gluon self-energy. This we 
do in order to improve the large Q^ behaviour of the coefficient function. 
Without this vertex correction L^ ' ~ ln^(Q^/m^) whereas it behaves like 
In (Q'^/m^) if this correction has been included. The coefficient function 
Hj^ ' is computed from the Bethe-Heitler process given by the reaction 



W + qiq)^qiq)+cic)+q'iq') (2.11) 



with q = u,d,s and q' = d, s. Like reaction (|2.1CI| ) process (|2.11|) shows up for 
the first time if the QCD corrections are calculated up to order a^. Contrary 



to ( |2.10| ), where the W-boson is only attached to the light quarks q and q', 



the vector boson in reaction ( 2.11| ) is now coupled to the charm quark and the 



light quark q'. Finally we have the W-boson-gluon fusion mechanism which 
contributes to H/^g . In order as (lowest order) this production mechanism 
is given by the process (see also ( |1.5| )) 



W + g^c{c) + q'{q') (2.12) 

with q' = d, s. In NLO we have to include all virtual corrections to reaction 
( |2.12| ) and to add the contributions due to the gluon bremsstrahlung process 



W + g^cic)+q'iq')+g. (2.13) 



Contrary to L/^ ' for which the exact order a^ expressions exist, see [IT? 



''k,q 

{k = 2) and ^ {k = 3) we do not have the exact form of H^f^ and H^f 
in the same order. This is in contrast to the heavy-flavour electroproduction 
where the deep-inelastic process only proceeds via the exchange of a virtual 
photon. In this case the exact order a^ contributions to the heavy-quark 
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coefficients H]^^ and if 2,'g ^^^ known and they can be found in [|T^ . For the 
charged- current process the calculation of these coefficient functions is still 
more difficult than the ones for the electromagnetic-current reaction as given 
above. This can be attributed to the mass of the light quarks q' in ( p.ll|) - 
( |2.13|) which is usually put to be equal to zero. Therefore the power of the 
coUinear singularities, occurring in the higher order radiative corrections to 
the charged-current reaction, is much higher than the one appearing in the 
electromagnetic-current process where also q' stands for the charm quark. 
This in particular will complicate the calculation of the many particle phase 
space integrals. The latter are even more complicated than those encountered 
in parton subprocesses where all masses of the particles in the final state are 
equal to zero. These type of phase space integrals have been calculated in 
p^ . However these complications can be avoided if one is only interested 



in the asymptotic form of the heavy-quark coefficient functions given for 
Q^ ^ m'^. This form can be derived using the renormalization group and 
mass-factorization techniques. These methods have been applied to obtain 
the heavy-quark coefficient functions in the limit Q^ 3> m^ for unpolarized 



[|T^ and polarized |]T8| electromagnetic-current deep-inelastic lepton-hadron 
scattering processes. Here in the FFS the light parton contribution to the 
structure function F2^{x, Q"^) (/ = e, /z) is given by 



F^^ix,Q')^F^^ix,Q 



=-i^3[-,n']C^!{z,^,') 




C(^'/^' 



(2.14) 



and the charm component reads 



mx,Q^m2)^F£(x,Q 
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2,9 P'™2' ,,2 


4 /-^th (iz 


\ n\ -^ ,,2 t77,S ^ V >"'' \ 

\z j '^V m2 /i2y 





(2.15) 

The notations above are the same as those given for the charged-current 
structure functions in ( p.l|) and ( p.7|) . Notice that the hght-parton coeffi- 
cient functions Cfc,j are the same in both reactions which means that they do 
not depend on the probe W or 7. This does not hold for the heavy-quark 
coefficient functions. The most important result of these calculations is that 
the F2^J^x, Q2) does not change when the exact heavy-quark coefficient func- 
tions are replaced by their asymptotic expressions provided Q^ > Ql^in- In 
the case of charm production Q^j„ = 20 (GeV/c)^ for x < 0.1. Further- 
more it turns out that the value obtained for Q^j„ is not altered while going 
form the LO to the NLO approximation. If we now assume that the same 
property also holds for the charged-current structure functions F^^{x,Q'^), 
Fl^^{x, Q"^) then the asymptotic forms of if^ can be used to make fairly good 
predictions of the order a1 contributions. These predictions will become very 
accurate for experiments carried out at HERA since the Q^- values are very 
large and the x-values are reasonable small i.e. 200 < Q^ < lO'' (GeV/c)^ 
and 0.006 < x < 0.5. Even for fixed-target experiments where the charac- 
teristic Q^-values are much smaller one can make reliable predictions unless 
one enters the threshold region where Q^ is too small [Q^ < 20 (GeV/c)^) 
and/or x is too large (x > 0.1). 

The asymptotic heavy-quark coefficient functions H^^ for the charged- 
current process are derived in Appendix A. It appears that they can be 
expressed into the asymptotic heavy-quark coefficient functions H^^ in ( p.l5| ) 
computed for the electromagnetic-current process (see Appendix D of ||T7|). 
Let us first first expand these coefficient functions in a power series like 



rV _ Y- / «M H-^'(0 




HL = Uiz H^'r (2.16) 



with V = J, W; k = 2,3 and i = q,g. From Appendix A it follows that the 
coefficients Hf, j up to order a^ in the limit Q^ ^ m? are related to each 
other as follows 
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Q2 rr^' 

9 ' 9 



rrW,S,il) I ^- "^- \ _ 1 rr7,S,{l) , 
^k,9 I ^'^2' ,,2 j - 2^'^ 



9 ' 9 



^'^'^l+^^'^Clf (^'^)'(2-17) 






i/, 



VK,S,(2) I 



k,g 






-h: 



7,S,(2) I 



2,9 



Q2 ^2 

m^' ^2 




4r^^"&f){^K"(.f)-ci-(.f)} 



(2.18) 



ttW,PS,{2)I Q^ ^ 

with 



2 ^•'^ 



7,PS,(2)/ ^Q^^rn^^ 



r^fcC; 



PS, (2) 
2,g 



>f) P^i^' 



62 = 1, 



es 



-1 



(2.20) 



where i7- 



7,(0 



2,i 



g, (yf) can be found in Appendix D of [|T7 



Further A^'^^^ 

stands for the one-loop heavy-quark operator matrix element defined in (A. 6) 
and Cfc^j (fc = 2, 3; i = g, (?) denote the light-parton coefficient functions which 
are presented up to order a\ in Appendix B of pSf. Here we want to empha- 



size that relations ( |2.17|) - (p.l9| ) only hold in the asymptotic limit Q^ ^ rn? . 
Further they are presented in the MS-scheme. One can check the validity 
of ( p.lTp using the exact order a^ contributions to the coefficient functions 
if^ and 7^2,9 computed in 0. Finally notice that the second order coeffi- 

cient functions Q / require the knowledge of the three-loop DGLAP splitting 
functions in order to make a complete next-to-next-to-leading order (NNLO) 
analysis of F^\ in ( |2.7| ). However these splitting functions are not known 
yet so that one cannot obtain the parton densities in the next-to-next-to- 
leading log (NNLL) approximation. Hence the order a^ contributions to the 
structure functions F^^ due to the expressions in ( |2.18| ), ( p.l9| ) have to be 
considered as an estimate of the exact QCD corrections beyond NLO. How- 
ever in practice it turns out that the contributions due to the coefficient 
functions are much more important than those coming from the higher order 
DGLAP splitting functions provided both are represented in the MS-scheme 
(see @). 
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3 0(a'^) contributions to the charged- current 
process due to charm production 

In this section we will discuss the QCD corrections to the charged-current 
structure functions F^dx, Q^, rrP') {k = 2, 3) measured at fixed isoscalar tar- 
get experiments : {v^ij)^) + N -^ fi~{fi~^) +'X') [0]-|Q and the collider ex- 
periments at HERA : (e"(e+) +p^ u^iue) +'X') fTl-IHl- The kinematical 



ranges in which the fixed-target experiments are carried out are given by 
0.01 < X < 1 and 2 < Q2 < 100 (GeV/c)^. For HERA they become 
0.006 < X < 0.5 and 200 < Q^ < 10^ (GeV/c)^. In particular we are in- 
terested in the order a^ contributions originating from the heavy-quark co- 
efficient functions Lf,'^ '^ (exact) and i^^^ , Hf,^ '^ (asymptotic) pre- 
sented in the previous section. Since we are dealing with approximations to 
the above contributions only we have to determine the kinematical range in 
which they are valid. First we have to investigate for which x- and Q^-values 
the order a^ corrections coming from H)^' (W- boson-gluon fusion pro- 
cess) dominate those which are due to Hj^^^ '^ (fiavour-excitation process). 
Then we have to determine Q^in so that for Q^ > (5^j„ the results obtained 
from the asymptotic and exact expressions of the heavy-quark coefficient 
function Hj^ g coincide. Because of the observations made in [|1^]-^0| for 
the one-photon exchange process we expect that the same values of Q^m 
also holds for the order a^ contributions to the charged-current interactions 
given by H^^ '",H,^^g'". 

In our plots we will use the notations that F^^ (/ = 0, 1, 2) receive contri- 
butions from the coefficients corrected up to order a^. This implies that for 
Ff, I (Born- approximation or LO) we have to adopt the leading log (LL) par- 
ton densities and the running coupling constant. In the case of F^j (NLO) 
we have to choose the next-to-leading log (NLL) parton densities and the 
running coupling constant. This will be also done for F^\: where we only 
have included the order a^ contributions due to the heavy-quark coefficient 
functions L^, ^ \ Hj^J^ '^ and Hj^^ '^ . All parton densities, coefficient 
functions and the running coupling constant are presented in the MS-scheme. 
Further we choose the factorization scale equal to the renormalization scale 
and set /i^ = Q^. For other scales and the scale dependence of F^c see 
[0, |]TB[. The structure functions in ( |2.7| ) are presented in the fixed-flavour 
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scheme (FFS) where the number of hght flavours is chosen to be nf = 3. 
Since the parton density set GRV94 [^| is the only one which is consistently 



determined in the FFS we will adopt this in our paper. In LO and NLO 
the QCD scale A„^., which has to be substituted in the running coupling 

constant, is given by A3 = 248 MeV and A3 = 232 MeV respectively. 
The remaining input parameters are the charm-quark mass and the Cabibbo 
angle which are chosen tohe m = 1.5 (GeV/c) and | sin^d = 0.221 |2^. For 
the dependence of F^c on m see 0. 

In order to determine the kinematical range in which the approximations 
are valid we will first plot the quantities F^^ and Ff^ which receive the 

order a^ contributions from if^ ' and Hf^' respectively. In Fig. 1 

we have presented F2C, F| ^ for an isoscalar target in neutrino reactions at 
Q^ = 10, 100 (GeV/c)^. Here one observes that the W-boson-gluon fusion 
process {F2^^) dominates the flavour-excitation process (-Fj'^c) when x < 0.07 
(Q2 = 10(GeV/c)2). At larger Q^ i.e. (q2 _ ^qO (GeV/c)^) x increases 
and the latter becomes equal to x = 0.2. The same plots are given for F3 ,, 
and F|c in Fig. 2. It turns out that for F^c the x-values, below which F^,, 
dominates F3 ^, are larger that the ones obtained for F2 c- In Fig. 2 they 
become x < 0.3 for both Q^ = 10 (GeV/c)^ and Q^ = 100 (GeV/c)2. From 
the above we conclude that already at rather low Q^-values the W-boson- 
gluon fusion mechanism completely dominates the order as corrections to 
-p2,c provided x < 0.1. This also holds for F^^c at even larger x- values i.e. 
X < 0.3. 

Now we have to determine at which value of Q^, denoted by Q^m? ^^e 
exact heavy-quark coefficient function Hj^^ '^ can be replaced by its asymp- 
totic expression without altering the structure functions F^ c- For that pur- 
pose we will study the quantity 



2\ k,c 



Fi:r'\x,Q^m^^ 



RkA^^Q ,rn) = —^j—^^i——- — (3.1) 



where the numerator and the denominator are computed using the asymp- 
totic and the exact expressions for H^^ respectively. Further we focus our 
attention on the range 0.006 < s < 0.1 which is explored by the fixed-target 
experiments and the HERA collider. The range 0.1 > a; < 1 falls beyond our 
scope because here the P^-boson-gluon fusion process is not dominant any- 
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more. In Fig. 3 we have plotted i?2,c in the range 10 < Q^ < 10^ (GeV/c)^ 
for X = 0.1, 0.01, 10~^, 10~^. The same is done for R^c in Fig- 4. Fig. 3 
reveals that for x = 0.1, R2,c > 0.9 when Q'^ > Qf^^^ = 50 (GeV/c)^. 
At decreasing x the value of (5^j„ will become less so that for x = 10~^ 
one obtains Q^j„ = 20 (GeV/c)^. A similar behaviour is shown by Fig. 4. 
Here we observe that for x = 0.1, R^^c > 0.9 when Q^j„ = 85 (GeV/c)^. 
Also in this case Q^m decreases if x gets lower. For instance for x = 10~^ 
we get Q'^in = 50 (GeV/c)^ . From the above one cannot conclude that 
the replacement of the exact by the asymptotic coefficient functions works 
better for F2 c than for F^c- This becomes clear when Fig. 5 is compared 
with Fig. 6 where now Rk^c ( p.l|) {k = 2, 3) is plotted as a function of x 
for Q2 = 10,20,50, 100 (GeV/c)^. If we choose Q^ = 10 (GeV/c)^ then in 
the range 0.01 < x < 0.1 we observe that 0.6 > R^c > 0.5 (Fig. 6). How- 
ever the behaviour of i?2,c is terrible. In the range 0.01 < x < 0.05 we find 
0.7 > -R2,c > 0.1. For x > 0.05, i?2,c becomes even negative and for very large 
X it switches sign again. This behaviour of the coefficient functions originates 
from mass factorization of the original partonic cross sections which destroyes 
the positive definiteness of the latter quantities. Further we also studied the 
behaviour of F^^, F^^ and Rk^c for antineutrino isoscalar reactions but the 
results are the same as found above. Summarizing our findings we can state 
that for X < 0.1 and Q^ > Qlnin = 50 (GeV/c)^ the approximation to the 
structure functions, given by the replacement of the exact by the asymptotic 
heavy-quark coefficient functions, works quite well. One can also choose 
smaller values for Q^m- However in this case the approximation will still 
work for F^c but it becomes bad for F2^c unless one chooses smaller values 
for X. For example for Q^j„ = 20 (GeV/c)^ one has to take x < 0.08. From 
now on we will assume that the same Q^m ^^^ ^ ^^so hold for the order a1 
contributions to the heavy-quark coefficient functions H^J^ , Hf,g This as- 
sumption is based on our findings in [|l^-[^ for the electromagnetic-current 
process where the above contributions are exactly known up to order a^. 
In the next figures we show ^2^ (/ = 0,1,2) for the neutrino isoscalar 

target reaction. Here we have plotted ^2,. (exact), Fj,, (exact) and F2J (ap- 
proximate) for g2 = 10 (GeV/c)2 (Fig. 7) and Q'^ = 100 (GeV/c)^ (Fig. 8). 
Fig. 7 reveals that the order a"^ contributions, which are positive, become 
very small. However here one has to bear in mind that for x > 0.03 and 
Qmin = 10 (GeV/c)^ the estimate of this contribution is not very rehable 
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in view of what we found in Fig. 5. At Q^ = 100 (GeV/c)^ (Fig- 8), where 
the approximation is very accurate, the order a^ corrections are larger and 
at a; = 0.01 they amount to 5% of F2J which is still quite small. These 
corrections in the case of ^3^ (see Figs. 9,10) become much larger than the 
ones observed for F2^c- Furthermore in the light of the discussion above for 
i?3^c (Fig. 6) they are more reliable even for Q^ = 10 (GeV/c)^ (Fig. 9). 
Like in the case of F2^c the largest corrections occur at the lowest x-value 
and they decrease when x tends to one. At x = 0.01 the order «^ cor- 
rections w.r.t. Fgc amount to 14% for Q^ = 10 (GeV/c)^ (Fig. 9) and 
20% for Q2 = 100'(GeV/c)2 (Fig. 10). Therefore at increasing Q^ the cor- 
rections get larger. Finally we want to call attention to the difference in 
behaviour between F2C and F^j, for / = 0, 1, 2. From Figs. 7-10 we infer that 

Fg > F« > Fi? whereas FJ? > F^? > F^^l 

We will now investigate how accurately one can determine the strange- 
quark density s{x, /x^) from the data obtained for the charged- current [Q] and 
electromagnetic-current |^ interactions. To that order we will compute the 
ratio in the three-flavour scheme which is defined by 



-i"(-«^'"^)-^SSr '^■^' 



with (see also (|1.6|)) 



JFf^^'^+FfJ^\ (3.3) 



where Fg^, F^^ (/ = u^u^jj) denote the light-parton contributions (( p.l| ), 
(|2.14|) ) and the charm component (( p.7|) , (p.l5|) ) of the structure functions 



respectively. For the computation of F!^^^ we will use the exact order a^ 
contributions to the heavy-quark coefficient functions as presented in [^ 



Further s*^°^(x, jj?) and s*^^^(x, jj?) stand for the LL and NLL parametrizations 
of the strange-quark density. Since the NNLL parametrization is unknown 
we have to put s*^^)(a;, /x^) = s'^^^(a;, /x^). The phenomenon that Rf"* f) can be 
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only due to the effect of tlie mass of tlie cliarm quark. For / > 1 one gets in 
addition deviations wliich originate from tlie higher order QCD corrections. 
We have plotted R^''^ (/ = 0, 1, 2) in the range 10"*^ < x < 1 for three different 
g2-values i.e. Q^ = 10, 100, 10^ {GeV/cf. In Fig. 11 (Q^ = 10 (GeV/c)^) we 
observe that R^^^ strongly deviates from one in the region 0.1 < x < 1, which 
is wholly due to thresholds effects of the charm quark in Fj^^ and F2C . 
This effect will vanish when Q^ gets larger (see Figs. 12,13). Notice that in the 
region x > 0.1 the size of the strange-quark density is very small. For x < 0.1 
R^^^ tends to one for Q^ > 10 (GeV/c)^. This picture does not change when 
one includes the order as corrections provided R^-^^ is computed at small Q^ 
(see Fig. 11). At large Q^ (see Figs. 12,13), R^^' differs apprecibly from one in 
particular at small x. The difference becomes even larger when one includes 
the order a^ contributions (see R^^^ in Figs. 11-13). This effect is wholly due 
to the heavy-quark coefficient functions if]^j {V = 7, PF), appearing in the 
structure functions of ( p.3| ) which grow as ln*((5^/m^) when Q^ ^ m^. Notice 
that from Fig. 5 one can see that the estimate of the order a^ contributions 
to Fl^ and Fl^ is quite reliable for Q^ > 10 (GeV/c)^ and x < 0.01. From 
Figs. 11-13 we can conclude that at small x, R^^^ and R^^^ become smaller 
than 0.5, which means that the actual strange-quark density is more than 
two times larger than the value obtained for A in ( ^.3] ). Hence one can 
conclude that the higher order QCD corrections bedevil the extraction of the 
strange-quark density in the region x < 0.1. Furthermore the effect of the 
mass of the charm quark is not negligible at small Q^ {Q"^ < 10 (GeV/c)^) 
in the large x-region {x > 0.1). 

The steep rise of the heavy-quark coefficient functions which grow like 
ln*(Q^/m^) for Q^ ^ m^ becomes even more conspicuous if we study charm 
production at HERA where one measures the structure functions F^''^{= 

Fl/) and F'^f{= F^^). The latter are plotted for Q^ = 10^ {GeV/cf 
in the region 0.006 < x < 0.5 in Figs. 14-17. In this kinematical region 
the approximations to the second order heavy-quark coefficient functions 
become excellent (see Figs. 3-6). This already happens at the lowest Q^- 
value i.e. Q^ = 200 (GeV/c)^ measured in ||16|. In the case of F2^c the exact 
order a^ corrections are not that spectacular and they amount to almost 
7% at x = 0.006 w.r.t. to Fgj (see Figs. 14,15). This is in contrast to 

F^^c (see Figs. 16,17) where the order a^ corrections become 32% w.r.t. F^J 
(see Fig. 17). These large corrections, occurring at small x and large Q^, 



vitiate the perturbation series in particular the one given for F^^c- Therefore 
the large logarithmic terms represented by hi^{Q'^ /m?) have to be resummed 
according to the VFNS approach (see ||TT|, |^). In the latter scheme the 
charm component of the structure functions takes the following form 



akix 



dz 





(3^4) 

where the parton densities and the light-quark coefficient functions Ck^i {k = 
2,3; i = q,g) are presented in the four-flavour scheme. The parton-density 
set now also contains the charm-quark density c{z,^'^) and the singlet com- 
bination reads 



S4(z, /i^) = £3(2:, /i^) + c{z, /i^) + c(z, /i^) (3.5) 

where T,^{z,p?) is defined in ( p.4|) . The structure function F^^ (= F^^ ^) 
can be derived by interchanging quarks and anti-quarks (see below ( p.9|) ). In 
order to compute (|3.4|) we need the charm-quark density which however is not 
present in GRV94 ||25| because the latter set is presented in a three-fiavour 
scheme. Therefore we have adopted GRV92 [^ with A4 = A4 = 

200 MeV which contains the charm-quark density. We first checked that using 
the densities in GRV92 for the expressions in the three-fiavour scheme one 
gets the same results for F2^c as obtained with GRV94. In the case of l-Psd, 
the latter set produces values which are 10% below those given by GRV92. As 
a comparison we have also plotted in Figs. 14-17 F^. j (VFNS) using Eq. ( |3.4| ) 

and GRV92. It turns out that the result for F^]} is larger in VFNS (^ than 
obtained by FFS (see Figs. 14,15). This observation was also made in [|T^]. 
The discrepancy becomes less when the order a1 corrections are included in 
the case of FFS (see F2J). This feature was also discovered for Fg^^^ in the 
electromagnetic-current process [EH]. However, for Fs^c the difference between 
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the two schemes is large (see Figs. 16,17), in particular in the small x-region. 
Here we observe that ^3^, in VFNS deviates appreciably from ^3^ as well 

as F^^ in FFS. This is no surprise because, as we have seen above, in FFS 
the large logarithms in the heavy-quark coefficient functions lead to larger 
corrections to F^c than to F2 c- Therefore the resummation of these terms 
will have a larger effect on F^^c than on F2_c- This is the explanation of the 
difference between the predictions of the FFS and VFNS approach for the 
HERA experiments as originally observed for F2J in [|l^]. Since we believe 



that the resummation of these large logarithmic terms has to be carried out 
for Q^ ^ m^, the results obtained by Eq. (|3.4|) (VFNS) are much more 
reliable than those obtained in FFS given by the expressions in Section |^. 

Summarizing our findings, we have studied the effect of the order a^ 
contributions from the heavy-quark coefficient functions to the charm com- 
ponent of the charged-current structure functions. Since for some of them 
the exact expressions do not exist we have to rely on approximations which 
are strictly speaking valid only for Q^ ^ m'^ (HERA-regime). However, from 
a comparison based on the exact and asymptotic order ag heavy-quark co- 
efficient functions corresponding to the W-boson-gluon fusion process, the 
predictions are also reliable for Q^ > 20 (GeV/c)^ and x < 0.08. At smaller 
X it is even sufficient to choose Q^ > 10 (GeV/c)^. 

The outcome of our analysis reveals that the order a^ corrections are small 
for F2^c but large for F3 c, in particular when Q^ is big. These large corrections 
are due to the logarithmic terms ln^{Q'^/m'^) dominating the heavy-quark co- 
efficient functions characteristic of the FFS approach. Therefore at large Q^ 
these terms have to be resummed in the context of VFNS, leading to a better 
prediction, in particular for F^^c- Finally we want to stress that the higher 
order QCD corrections bedevil the extraction of the strange-quark density 
from the charged- and electromagnetic-current total structure functions, in 
particular at small x {x < 0.08). 



20 



Appendix A 

In this section we will derive the relations between the heavy-quark coeffi- 
cient functions if^ and 112^ {k = 2,3; i = q,g) obtained for the charged- 
current and the electromagnetic-current processes respectively. Using mass- 
factorization techniques one has shown in [0-O that there exists a relation 



between the asymptotic (Q^ ^ rn'^) heavy-quark coefficient functions H^^ 
and the light-parton coefficient functions Ck^i- Hence one can derive that 
up to order a"^ the heavy-quark coefficient functions corresponding to the 
reactions 



i + V^qi + q2, (A.l) 

i + V ^qi + q2+J, (A.2) 



m 
in 12 



the asymptotic limit Q^ ^ m^ can be written as (see Eqs. (2.34), (2.35) 



k.i I 2 ' 2 ' 9 ' 9 




For convenience we have suppressed the Bjorken scaling variable z - 
and ® denotes the convolution symbol defined by 

{f®g){z) = I dzi I dz2 6{z - ziZ2)f\zi)g{z2) . (A.4) 



JO 



The light partons involved in reactions ( [A.l ) and (A. 2) are given by 



I 

and i {i,j = q,g). For V = 'y both gi and q2 stand for the heavy quarks 
whereas for V = W one of them is represented by a light quark. The masses 
of gi and g2 are given by mi and m2 respectively. The light-parton coefficient 
functions Ck^i (^ = q, g) are obtained from the massless-parton subprocesses 

^ + V^Jl+J2 + ...+Jk (A.5) 
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which are independent of the nature of the probe V. Further when the 
quarks gi and g2 in ( A. 3 ) are massless the coefficient functions satisfy the 



relation Cf ^^ = Cf g^- The hght parton coefficient functions in ( |A.3| ) are 
calculated up to order a'j. in [|23|. Expression ( [A.3|) requires the knowledge 
of the renormalized heavy-quark operator matrix elements (OME's) 



A, 



Qli I 



mr 



/i- 



<i\o,M^> 



(A.6) 



and the heavy-quark loop contributions to the renormalized light-parton 
OME's given by 



with 




(A.7) 



4?U^ =M(1-)- 



(A.8) 



The above OME's are calculated up to order a^ in [O and the renormalized 



expressions are presented in the MS-scheme in Appendix B of |^0[- Further- 
more they satisfy the conditions when mf = 



v4,,,(0) = 0, A,i{z,0) = 6,i6(l-z) 



(A.9) 



The local operators Oj{x) appear in the light-cone expansion of the product of 
the two electroweak currents which, after having taken the Fourier transform, 



show up in the expressions for the structure functions Fk{x,Q ) (see p8| 



In the case oi V = j Eq. ( [A .31 ) leads in the first and second order to the 
following relations {mi = m2 = rn^qi = q2 = c) 



Q 



^ m^' 



f77,(l) 

2,9 2' 2 




(A. 10) 



'^ \m?^ fj? J 



4?g).4j'(^; 



+c, 



S,(2) 
2,9 



'q_ 

.^^ 



^NS,(1) /Q^ 



•-2,9 



,/^ 



(A.ll) 
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where we have expanded all quantities up to order (a;^/47r)^. Notice that in 
the case of K = 7, ifg^ = 0. Furthermore Ci^c = Ci,q when m -^ 0. 

For V = W the above relations become {mi = m, m2 = 0; gi = c, q2 = 

d, s) 



For the heavy-quark coefficient functions H^^ we have to use the following 
relations 



Cr. = -Cll Cl = 0, ClJ = (A.16) 

which follow from charge conjugation invariance of the strong interactions. 
Hence we can derive from Eq. ( |A.3| ) 



<1i.'^)4<'(^)' <--' 



<'1S'^) 4-s'(^) 4-ni^) «^-"&) ■ '"^* 
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m 



rTW,PS,i2) ( Q^ ■■- 



1 .PS,(2)[^ 



(A.19) 



Similar relations can be derived when mi = 0, 1712 = m and gi = d, s; q2 = c. 
Using Eqs. ( |A.10|) -( [A.12| ) one can now express if^ into H'^i- The results are 
presented in (pl7|)-(pA9D. 
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Figure Captions 



Fig. 1. The order a^ contributions to F2^ ( |2.7| ) due to flavour excitation : 
F2^c ( |1-4| ) and ly-boson-gluon fusion : Ff^, (pT5|); 



solid line: F|_, at Q^ = iq (GeV/c)^, dotted line: ^2% at Q^ = iqo 
(GeV/c)2, dashed line: F|_, at Q^ = iq (GeV/c)^, dashed-dotted line: 
Fl^ at Q2 = 100 (GeV/c)2'. 

Fig. 2. Same as Fig. 1 but now for Fi^^ ( |2.7| ) plotted in absolute value. 



Fig. 3. -R2,c ( |3.1| ) plotted as a function of Q^ at flxed x; 

a; = 10~^ (dotted line), x = 10~^ (dashed-dotted line), x = 10~^ 
(dashed line) and x = 10~^ (solid line). 



Fig. 4. Same as Fig. 3 but now for R 



3,c 



Fig. 5. i?2,c (|3.1| ) plotted as a function of x at flxed Q^; 

Q2 =' 10 (GeV/c)2 (dotted line), Q^ = 20 (GeV/c)^ (dashed-dotted 
line), g2 = 50 (GeV/c)^ (dashed line), Q^ = 100 (GeV/c)^ (sohd line). 



Fig. 6. Same as Fig. 5 but now for R 



3,c 



Fig. 7. The order a^ corrected structure function F2^ ( p. 7] ), denoted by 
^2% as a function of x at Q'^ = 10 (GeV/c)^; 
dotted line : F2 1 , dashed line : F^ I , solid line : F2 1 . 



Fig. 8. Same as Fig. 7 but now at Q'^ = 100 (GeV/c)^. 

Fig. 9. Same as Fig. 7 but now for F^^ (|3) at Q^ = 10 (GeV/c)^ 

Fig. 10. Same as Fig. 7 but now for F^^ (|3) at Q^ = 100 (GeV/c)^. 
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Fig. 11. The order a[ corrected ratio Rg (|3.2|), denoted by R^^\ as a func- 
tion of X at Q2 = 10 (GeV/c)2; 
solid line : i?(°), dashed line : i?W, dotted hne : R'^\ 



Fig. 12 Same as Fig. 11 but now at Q^ = 100 {GeV/cf. 
Fig. 13 Same as Fig. 11 but now at Q^ = 10^ (GeV/c)^. 

Fig. 14 The order a^ corrected structure function ^2,,^ ( = F2^^) (|2.7|) , de- 
noted by F2 J, as a function of x at Q^ = 10^ (GeV/c)^; 
dotted line : F^^J (FFS), dashed line : fQ (FFS), sohd line : F^J 

(FFS). As a comparison we have also shown dashed-dotted line : ^2^, 
(VFNS). 



Fig. 15 Same as Fig. 14 but now for Fl/ ( = Fs^f ) (g 
Fig. 16 Same as Fig. 14 but now for F3/ ( = F2';f ) (| 
Fig. 17 Same as Fig. 14 but now for F^/ ( = ^2^^ ) (|2 
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